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Abstract 

A set of closed form equations for the preliminary evaluation and design of a heat- 
pipe-cooled leading edge is presented. The set of equations can provide a leading-edge 
designer with a quick evaluation of the feasibility of using heat-pipe cooling. The heat 
pipes can be embedded in a metallic or composite structure. The maximum heat flux, total 
integrated heat load, and thermal properties of the structure and heat-pipe container are 
required input. The heat-pipe operating temperature, maximum surface temperature, heat- 
pipe length, and heat pipe-spacing can be estimated. Results using the design equations 
compared well with those from a 3-D finite element analysis for both a large and small 
radius leading edge. 


Nomenclature 


English 

A 

B 

h 

k 

L 

P 

q 

q” 

r 

R 

s 

t 

T 

w 

X 

y 

z 


area, in 2 

parameter utilized for grid transformation, in. 

heat transfer coefficient for contact resistance, Bhi/ln--ft 2 -°F 

thermal conductivity, Btu/hr-ft-°F 

heat-pipe length, in. 

equally spaced coordinates 

heat flux, Btu/s 

heat flux per unit area, Btu/ft 2 -s 
leading-edge radius, in. 
thermal resistance, hr-ft 2 -°F/Btu 
heat-pipe length, in. 
thickness, in. 
temperature, °F 
width of heat pipe, in. 
half heat-pipe spacing, in. 

transformed coordinate for grid transformation, in. 
point about which grid clustering occurs, in. 


Greek 

£ emittance 

o Stefan Boltzmann constant 

2 summation 

x stretching parameter for grid transformation 



Subscripts and superscripts 
amb ambient 

anal obtained from 1-D closed-form design equations 

avg average 

c coating 

cr contact resistance 

FEA finite element analysis 

hp heat pipe 

L lower 

max maximum 

s,hp difference between structure and heat pipe temperature 
s,t structure property, through-the- thickness direction 

s,p structure property, in plane direction 

stag stagnation 

surf surface 

tot total 

U upper 

w heat-pipe wall 


Introduction 

Stagnation regions, such as wing and tail leading edges and nose caps, are critical 
design areas of hypersonic aerospace vehicles because of the hostile thermal environment 
those regions experience during flight. As a hypersonic vehicle travels through the earth's 
atmosphere, the high local heating and aerodynamic forces cause very high temperatures, 
severe tlienual gradients, and high thenual stresses. Analytical studies, laboratory, and 
wind tunnel tests indicate that a solution to the thenual- structural problems associated with 
stagnation regions of hypersonic aerospace vehicles might be obtained by the use of heat 
pipes to cool these regions. 

In the early 1970's, several feasibility studies were performed to assess the application 
of heat pipes for cooling leading edges and nose caps of hypersonic vehicles. 1 NASA 
Langley Research Center (LaRC), through a contractual study, analytically verified the 
viability of heat pipes for cooling stagnation regions of hypersonic vehicles. 1 In 1972, 
McDonnell Douglas Astronautics Co. (MDAC) compared four space shuttle wing leading- 
edge concepts: a passive carbon-carbon concept, a passive coated-columbium concept, an 
ablative concept, and a liquid-metal/ superalloy heat-pipe-cooled concept. 2 The heat-pipe- 
cooled concept was determined to be a feasible and durable design concept, but was 
slightly heavier than the other candidate concepts, hi 1973, MDAC fabricated a half-scale 
shuttle-type heat-pipe-cooled leading edge to verify feasibility of the concept. 4 This model 
was tested by a series of radiant heating and aerothermal tests at NASA LaRC from 1977 to 
1978 to verify heat-pipe transient, startup, and steady-state performance. 6 8 hi 1979, 
MDAC received a follow-on contract to optimize a heat-pipe-cooled whig leading edge for a 
single- stage-to- orbit vehicle. Results of the follow-on study indicated that the mass of a 
shuttle-type heat-pipe-cooled leading edge could be reduced by over 40% by use of a more 
efficient structural design. 9 hi 1986 MDAC received a contract to design and fabricate a 
sodium/superalloy heat-pipe-cooled leading edge component for an advanced shuttle-type 
vehicle. 10 This advanced shuttle- type heat pipe was 6-ft long and was tested at MDAC by 
radiant heating and at Los Alamos National Laboratory (LANL) by induction heating. 11 12 

Preliminary design studies at NASA LaRC indicate that a refractory- 
composite/refractory-metal heat-pipe-cooled leading edge can reduce the leading-edge mass 
by over 50% compared to an actively cooled leading edge, can completely eliminate the 
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need for active cooling, and has the potential to provide failsafe and redundant features. 13 
Recent work to develop this novel refractory-composite/refractory-metal heat-pipe-cooled 
leading edge for hypersonic vehicles combines advanced high-temperature materials, 
coatings, and fabrication techniques with an innovative thermal- structural design. Testing 
of a component at NASA LaRC with three straight molybdenum-rhenium (Mo-Re) heat 
pipes embedded in carbon/carbon (C/C) has demonstrated the feasibility of operating heat 
pipes embedded in C/C. 14 15 

When confronted with a leading-edge design for hypersonic vehicles, the options are 
passive, heat-pipe cooled, or actively cooled. The upper use limit for passive leading edges 
may be determined by evaluating the material properties hi light of the thermal and 
mechanical loads. If passive leading edges cannot survive the environmental conditions, 
heat-pipe cooled or actively cooled leading edges will be required. Though heat pipes are 
often a viable and light weight option, the analysis required to detennine if heat pipes are 
feasible for a particular application can be extensive and may thus preclude their use. It is 
thus beneficial to have a simple set of closed fonn equations that can be used to detennine if 
the heat-pipe option is feasible. Having a simple analysis available may prevent the 
unnecessary use of active cooling when heat-pipes may provide a cheaper and lighter 
weight alternative and may prevent the unnecessary use of complex, 3D finite element 
analysis techniques to answer the question of initial feasibility, thus saving substantial 
analysis time. 

The purpose of this paper is to present a set of simple, closed-fonn design equations 
that can be used to determine the feasibility of using a heat-pipe -cooled leading edge. The 
design equations presented here are only for thennal design, and do not include any stress 
analysis. Temperatures obtained from the design equations are compared to a 3-D finite 
element analysis for both a large and small leading-edge radius. Though some restrictions 
apply to the use of the equations, they appear to be a useful tool for a preliminary look at 
the feasibility of heat-pipe-cooled leading edges. If the preliminary design equations 
indicate a feasible design, a more detailed analysis should follow. 

Description of Heat-Pipe-Cooled Leading-Edge 

A brief description of how heat pipes operate and are utilized for leading-edge cooling 
is first presented, followed by a brief description of the heat-pipe-cooled leading-edge for 
which the equations were developed. 

Leading-Edge Heat-Pipe Operation 

Heat pipes transfer heat nearly isothennally by the evaporation and condensation of a 
working fluid, as illustrated in Figure 1. The heat is absorbed within the heat pipe by 
evaporation of the working fluid. The evaporation results in a slight internal pressure 
differential that causes the vapor to flow from the evaporator region to the condenser 
region, where it condenses and gives up heat. The cycle is completed with the return flow 
of the liquid condensate to the evaporator region by the capillary action of a wick. 

Heat pipes provide cooling of stagnation regions by transferring heat nearly 
isothennally to locations aft of the stagnation region, thus raising the temperature aft of the 
stagnation region above the expected radiation equilibrium temperature. When applied to 
leading-edge cooling, heat pipes operate by accepting heat at a high rate over a small area 
near the stagnation region and radiating it at a lower rate over a larger surface area, as 
shown hi Figure 2. The use of heat pipes results in a nearly isothermal leading-edge 
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surface, thus reducing the temperatures in the stagnation region and raising the 
temperatures of both the upper and lower aft surfaces. 
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Figure 1: Schematic diagram of the operation of a heat pipe showing the heat-pipe 

container, working fluid, and wick. 



Figure 2: Schematic diagram of a heat-pipe-cooled leading edge showing regions of net 

heat input (evaporator) and net heat output (condenser). 

Refractory Composite Heat-Pipe-Cooled Wing-Leading-Edge 

The refractory composite heat-pipe-cooled wing leading edge for which the design 
equations were developed is illustrated in Figure 3. The heat pipes are oriented normal to 
the leading edge and have a “D-shaped” cross section, with the flat paid of the “D” forming 
the wing-leading-edge outer surface. As shown in Figure 3, the leading edge contains “J- 
tube” heat pipes, with a “J-tube” heat pipe being a heat pipe with a long leg on one side of 
the nose region, and a short leg on the other side of the nose region. An alternating “J- 
tube” configuration was selected here to minimize heat-pipe spacing in the nose region 
where heating is the highest, provide a greater heat-pipe spacing on the upper and lower 
surfaces where heating is lower, and, at the same time, minimize mass. The refractory 
composite structure sustains most of the mechanical structural loads and also offers ablative 
protection in the event of a heat-pipe failure. 

The maximum operating temperature capability of coated refractory-composite 
materials for the primary structure of the leading edge is high (~3000°F) relative to 
refractory metals, which are typically limited to approximately 2400°F. The potentially 
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higher operating temperature increases the radiation heat-rejection efficiency of the heat- 
pipe-cooled leading edge and permits reductions in the mass of the leading edge for a given 
leading-edge radius. In addition, the higher operatmg temperature increases the total heat 
load that can be accommodated passively by the heat pipe (i.e., no forced convective 
cooling required). For many trajectories, the high operatmg temperatures help eliminate the 
need for active cooling during both ascent and descent, thus eliminating the need for 
carrying additional hydrogen fuel (coolant) mto orbit. Since many hypersonic vehicles 
return unpowered for landing, the additional hydrogen fuel needed for cooling during 
descent would result in a mass penalty. 



Figure 3: Schematic drawhig of a hypersonic vehicle with a diagram of a heat-pipe-cooled 
wing leading edge. 


Design Equations 

The design of a heat-pipe-cooled leading edge is very complex due to the numerous 
variables involved. However, a simple set of closed form equations is presented here that 
can be used to determine if a heat-pipe-cooled leading edge is feasible with various material 
combinations. The equations presented here were developed to model the heat-pipe-cooled 
leading edge shown in Figure 3, but can be generalized for many other potential designs. 

Applied heat flux 



Figure 4: Schematic drawing of three heat pipes embedded in a structure. 

Figure 4 shows a schematic cross-section diagram of three heat pipes embedded in a 
structural material. The heat pipes shown in the figure have a rectangular cross section, but 
other cross sections could be considered. The leading edge is subjected to aerodynamic 
heating on the outer surface. At the stagnation line - the location of maximum heating - the 
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heating rate is by denoted q” stag . A coating of thickness t c is placed on the outer surface of 
the structure. The thickness of the structure between the outer surface and the heat pipes is 
t s and the heat-pipe container wall thickness is t w . The distance between heat pipes is 2x, 
and the width of the heat pipe is w. Contact resistance between the structure and the heat 
pipe is also shown in the figure. The contact resistance on the other surfaces of the heat 
pipe is of much less concern and is thus neglected in this calculation. 


The first step is to detennine the temperature drop, AT stag , through the structure and 
heat-pipe container at the stagnation line. This will help detennine the maximum 
temperature of the leading edge, which will occur on the outer surface at the stagnation line 
midway between heat pipes. To detennine the maximum temperature drop through the 
structure and heat-pipe container at the stagnation line, the following thennal resistances 
should be considered: through-the-thickness of the structure (from the outer surface to the 
heat pipe), in the plane of the structure (from midway between heat pipes to the heat pipe), 
and the contact resistance. If a coating is used on the outer surface, its thennal resistance 
(both in-plane and through-the-thickness) should be included. Two conduction paths are 
shown in Figure 5 for the heat conducted from midway between heat pipes on the outer 
surface to the heat pipe. As shown hi Figure 5, the heat must be conducted through the 
coating and structure in the through-the-thickness direction, and through either the coating 
or structure in the in-plane direction. 



Structure 


Figure 5: Schematic drawing of leading-edge cross section with heat pipe showing two 

potential paths for heat to be conducted from midway between heat pipes to a 
heat pipe. 


Heated surface 



Through-the-thickness 


In plane 


Figure 6: 


Schematic drawing of the series/parallel resistance network for heat to be 
conducted from midway between heat pipes on the outer surface to a heat pipe. 
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The drawing shown hi Figure 6 is an attempt to approximate the 2-D geometry of 
Figure 5 with a 1-D thermal resistance network. The thermal resistance is for the heat 
conduction midway between heat pipes on the outer surface to the heat pipe. The through- 
the-thickness resistance in the structure is shown in Figure 6 prior to the in-plane 
resistances, but could be placed after the in-plane resistances with the same result. Other 
resistance networks could also be used, but care must be exercised due to inconsistent 
areas. The thermal resistance through the heat-pipe container is neglected since it is small 
relative to the other tenns and since the heat transfer area is not consistent with the other 
tenus. Contact resistance is not shown hi Figure 6 due to its unknown value and area. 

Knowing the stagnation heat flux, q” stag , the dimensions, and the thennal 
conductivities, the temperature drop from a point midway between heat pipes on the outer 
surface to the heat pipe, AT stag , can be determined from 


q” st ag = AT stag /IR 


AT 


stag 


stag 


_^ + JiL + 
k c k s ,t L 


k sp t s /(t s +t c ) | k c t c /(t s + t c ) ~ 


(la) 

(lb) 


The thermal resistance given in eq. (lb) is for the geometry shown hi Figure 6, which was 
an attempt to approximate the 2-D geometry in Figure 5 with a 1-D thermal resistance 
approach. The first two tenus represent the through-the-thickness series resistance through 
the coating and structure. The third term represents the in-plane parallel resistance through 
the coating and stmcture. The thickness ratios, t s /(t s + t c ) and t c /(t s + t c ), represent the cross 
sectional area (assuming a unit depth) for the heat conduction through each layer. 
Rearranging eq. ( lb), the temperature drop can be obtained from 


AT 


stag 


= q 


stag 



x(t s +t c ) 1 

— + — + 

k c k s ,t 

Apt.+Mj 


( 2 ) 


The maximum value of the stagnation heat flux is used and the transient nature of the 
heating is not taken into account. Though this is a conservative approach, the thennal 
response of the leading edge will often be rapid enough that a steady state approximation at 
the time of maximum heating will provide relatively accurate temperatures. 

The next step is to determine the average surface temperature based on the estimated 
heat-pipe length. To do this, one must know the heat-flux distribution and estimate the 
chordwise length of the heat pipes on both the upper and lower surface. (The heat pipes 
will normally be oriented perpendicular to the leading edge, referred to as the chordwise 
direction, but could be oriented in the flow direction. For a swept leading edge, orienting 
the heat pipes perpendicular to the leading edge results in easier fabrication and lower axial 
heat-pipe acceleration loads.) From the heat flux distribution, the integrated heat flux, q tot , 
can be obtained for the entire chordlength, both upper and lower surface for a spanwise 
(parallel to the leading edge) unit width. The average outside surface temperature, T surf , can 
then be estimated from 


= e o A (T surf 4 - T^ 4 ) 


(3) 
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Thermal radiation to the leading edge from the ambient is included in eq. (3), but can 
usually be neglected. Reducing the upper or lower length of the heat pipes, L u or L, 
respectively, will raise the average surface temperature. The area, A, is based on a 1 -in- 
wide strip the total length of the heat pipe (L u + L,). It is important that the heat pipes 
extend past the stagnation region into the region where the maximum material reuse 
temperature is above the radiation equilibrium temperature. For sharp leading edges with 
small angles of attack, the heat flux drops off very rapidly, and should not be a problem. 
However, for very blunt leading edges with high angles of attack, high heat fluxes will 
extend a significant distance from the stagnation line. 

The thh'd step is to estimate the internal heat-pipe temperature. It is assumed that the 
heat pipe is at uniform temperature and that the heat radiated from the surface is also 
uniform. First, the heat flux out of the heat pipes is calculated assuming that the heat flux 
radiated from the leading-edge outer surface must first be conducted through the heat-pipe 
width, w. The distance between heat pipes is 2x and, thus, for every spanwise unit width 
of leading edge, the heat flux must be conducted through a heat-pipe of width w and is 
radiated from the outer surface over a width of w + 2x. Therefore, for each 1 in. unit 
width of leading edge, the heat is conducted to the outer surface through a width of 

(1 in.) w/(w + 2x) 

The average heat flux conducted through the wall of the heat-pipe container is then 


(L u + Lj) (1 in.) [w/(w + 2x)] 


Knowing the average heat flux, the temperature drop through the structure and heat-pipe 
container, AT shp , can be obtained from 



(5) 



Figure 7: Schematic drawing of a leading-edge cross section with a heat pipe showing the 

four thermal resistance components (coating, structure, contact resistance, and 
heat-pipe container) hi eq. (5). 


Note that in eq. (5) there is no in-plane thennal resistance. Contact resistance between the 
heat pipe and structure, if known, can be included in eq. (5). If values of contact resistance 
can be estimated, it may also be possible to bound the problem. However, it should be 
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emphasized that thermal contact resistance is extremely dependent on geometry, pressure, 
material, and temperature. 

The heat-pipe operating temperature is then obtained from 

Th P = T surf - AT shp (6) 

The final value to obtain is the maximum leading-edge temperature, which will occur 
midway between heat pipes at the stagnation line. Tins temperature is obtained from 

T max = T hp + AT stag (7) 

The important parameters for a heat-pipe-cooled leading edge have now been 
obtained: the maximum surface temperature, T niax , and the heat-pipe operating temperature, 
T hp . A comparison of the calculated maximum surface temperature with the reuse 
temperature of the coating and structural materials will determine if they are feasible for this 
application. The heat-pipe operating temperature will help establish what container material 
and working fluid to use. Several iterations may be required to obtain a design with 
acceptable temperatures using the corresponding material properties. The dimensions used 
in the design, i.e., the length of the heat pipes on the upper and lower surfaces, the spacing 
between heat pipes, the width of the heat pipes, and all the thicknesses, can be modified to 
obtain alternate designs. In addition, different materials with different thennal properties 
can be evaluated. The goal is to obtain a design that results in temperatures within the reuse 
limits of available materials while utilizing dimensions that can be fabricated. 

Comparison of Design and Finite Element Analysis Results 

The use of the developed design equations are now illustrated for a blunt leading 
edge, and results for both a blunt and sharp leading edge are presented and compared with 
results from a 3-D finite element analysis. 



Figure 8: Typical heat flux distribution for Space Shuttle Orbiter wing leading edge. 

As an example of the above procedure for a blunt leading edge, consider the Space 
Shuttle Orbiter wing leading edge with a stagnation heat flux of q’’ = 83 Btu/ft 2 -s and the 

chordwise (normal to the leading edge) heat flux distribution on the leading edge as shown 
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in Figure 8. Baseline (test case 1) dimensions and thennal conductivities are given in Table 
1. An example of a spreadsheet used for the calculations is shown in Appendix A. 

The temperature drop from the location at the outer stagnation line midway between 
heat pipes is given by eq. (2) as 


AT stag = q" 


stag 


— + — + 

x(t s +t c ) 

k c k s ,t 

>s,pt s +k c tcj 


( 2 ) 


AT 


stag 


83- 


Btu 


0.0 lin. 


ft 2 


23.587- 


Btu 
hr - ft - °F 


■ + ■ 


12 . 8 - 


0.04 in. 

Btu 
In- - ft - °F 


+ 


0.7in.(0.04in. + O.Olin.) 


^24.408- BtU 




hr - ft - °F 


(0.04 in.) + 


f 


23.587- 


Btu 


(2) 




hr - ft - °F y 


(O.Olin.) 


AT stag = 807°F 


Table 1: Baseline Variables for Comparison of Results 


t s 

= 

0.04 in. 

t w 

= 

0.01 in. 

t c 

= 

0.01 in. 

X 

= 

0.7 in. 

L u 

= 

24 in. 


= 

18 in. 

r 

= 

9 in. 

k s,t 

= 

12.8 Btu/hr-ft-°F 

k s,p 

= 

24.408 Btu/hr-ft-°F 

k W 

= 

39.485 Btu/hr-ft-°F 

k c 

= 

23.587 Btu/hr-ft-°F 

8 

= 

0.8 


The next step is to estimate the average leading-edge surface temperature. The length 
on the upper surface is 24 in. and on the lower surface is 18 in. The total integrated heat 
load for a spanwise width of 1 in. and a chordwise length of 42 in. is 1 1.6 Btu/s. Since all 
of the heat must be reradiated, and neglecting radiation from the ambient, the average 
surface temperature can be calculated from 


11.62 Btu/s = 8 o A T surf 4 = (0.8) (0.1714 x 10 8 Btu/hr-ft 2 -°R 4 ) (1 in.) (42 in.) T s 4 (3) 


T 


surf 


= 2738°F 


The heat-pipe container width is 0.6 in. with a heat-pipe spacing (half the distance 
between heat pipes) of 0.7 in. Thus for every 1 in. of leading-edge width (spanwise 
direction), 1 1.62 Btu/s must be conducted to the outer surface through a width of 


1 in. [w/(w+2x)] = 0.3 in. 
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The average heat flux conducted through the heat pipe to the outer surface is then 


11.62Btu/s 
avg ~~ (42in.)(0.3in.) 


132.8 Btu/ft 2 -s 


( 4 ) 


The thennal resistance from the coating, structure, and heat-pipe container are used to 
obtain the temperature drop through the leading edge. There is no in-plane thennal 
resistance included here since this is the region directly over a heat pipe. 


132.8 


Btu 

ft’s 


AT 


s,hp 


0.04in. 


12 . 8 - 


Btu 


+ - 


0.0 lin. 


hr - ft - °F 


39.485- 


Btu 


+ - 


0.0 lin. 


(5) 


hr - ft - °F 


23.587- 


Btu 


hr - ft - °F 


AT s , hp =151°F 


Table 2: Comparison of Blunt Leading-Edge Temperatures 


- 

r 

lip, anal’ 

T 

1 hpn-EA’ 

dp 

Difference, 

°F 

T 

max, anal’ 

°b 

T 

x max, FEA’ 

T 

Difference 

°F 

Test case 1 

q tot = 11.6 Btu/s 
t s = 0.04 in. 
x = 0.7 in. 

L u = 24 in. 

2587 

2556 

29 

3394 

3243 

150 

Test case 2 

q tot = 12.09 Btu/s 
t, = 0.04 in. 
x = 0.05 in. 

L u = 36 in. 

2530 

2420 

110 

2670 

2505 

165 

Test case 3 

q tot = 1 1.6 Btu/s 
t s = 0.25 in. 
x = 0.7 in. 

L u = 24 in. 

1933 

2525 

592 

3145 

3074 

71 

Test case 4 

q tot = 11.6 Btu/s 
t s = 0.25 in. 
x = 0.05 in. 

L u = 24 in. 

2455 

2537 

82 

3003 

2821 

182 


With an average outer surface temperature of 2738°F and a temperature drop through 
the structure of 151°F, the heat-pipe temperature is given as 


T hp = T surf - AT shp = 2738°F - 151°F = 2587°F (6) 
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The maximum leading-edge temperature is the sum of the heat-pipe operating temperature 
and the temperature drop from the location midway between heat pipes to the heat pipe, 
given by 


= T hp + AT stag = 2587°F + 807°F = 3394°F (7) 

A comparison of the heat-pipe temperature and maximum surface temperature are 
summarized in Table 2 for the 1-D design equations and a full 3-D finite element analysis 
(FEA) for a blunt leading edge such as on the Space Shuttle Orbiter. A discussion of the 
FEA is presented in Appendix B. Though the FEA has a nonlinear property capability, the 
constant properties in Table 1 were used in the FEA to provide a true comparison with the 
design equations. In each case, the variables that are different from those in Table 1 are 
listed in the table (all other variables are the same as in Table 1). For all cases in Table 2, 
q stag = 83 Btu/ft 2 -s and r = 9 in. 


Table 3: Comparison of Sharp Leading-Edge Temperatures 


T 

hp, anal’ 

°F 

T 

1 hnJFEA’ 

ip 

Difference, 

°F 

T 

max, anal’ 

T 

max, FEA’ 

°F 

Difference, 

°F 

Test case 5 2968 

q tot = 19.142 Btu/s 
t s = 0.04 in. 
x = 0.05 in. 

L u = = 24 in. 

3036 

68 

4230 

4077 

153 

Test case 6 2791 

q tot = 22.906 Btu/s 
t s = 0.04 in. 
x = 0.05 in. 

L u = Lj = 36 in. 

2851 

60 

3962 

3926 

127 

Test case 7 2826 

q tot = 19.142 Btu/s 
t s = 0.04 in. 
x = 0.7 in. 

L u = Lj = 24 in. 

2985 

159 

10,121 

6064 

4057 


Test cases 1-4 are for a blunt leading edge (r = 9 in.) with a relatively low heat flux 
and a large angle of attack (see Figure 8). Test case 2 has a much smaller half heat-pipe 
spacing than in test case 1 (0.05 in. vs. 0.7 in. in test case 1) and a longer upper surface 
heat-pipe length. The longer upper surface heat pipe results in a slightly larger integrated 
heat load. In test case 3, the thickness of the structure beneath the coating is increased to 
0.25 in. and the half spacing between heat pipes is 0.7 in. Both of these dimensions are 
relatively large and result in a heat-pipe temperature that is quite low. The combination of a 
thick structure above the heat pipe and a relatively large distance between heat pipes results 
in a large thermal resistance, and thus a large temperature difference, between the outer 
surface and the heat pipe. The larger dimensions also result in the 1-D approximation being 
less accurate, hi test case 4, the structural thickness is still large, but the heat pipes are 
spaced much closer, and the heat-pipe temperature from the design equations is much 
closer to that from the FEA. 
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Three cases are presented for a sharp leading edge (r = 0.5 in.) with a higher heat flux 
(q stag = 750 Btu/ft 2 -s ) in Table 3. The angle of attack for test cases 5-7 is near zero. The 
design results compare well with the FEA except for large heat-pipe spacing (test case 7), 
where the maximum temperatures obtahied by the two methods are very different. This is 
due to the fact that as the heat-pipe spacing increases, the problem becomes more three 
dimensional, and the design equations become less accurate. However, both the design 
equations and the FEA indicate that the design with a large heat-pipe spacing is not feasible. 

Discussion 

Since a heat pipe redistributes thermal energy instead of removing it as in active 
cooling, the total energy balance is extremely important. For this reason, sharp leading 
edges are much more conducive to heat-pipe cooling than blunt leading edges. A blunt 
leading edge, though it will have a lower stagnation heat flux than a sharp leading edge 
under the same flow conditions, may have a higher integrated heat load. As a result, the 
surface area required to radiate the energy away may be larger, i.e. longer heat pipes are 
required. 

Low angles of attack are more conducive to heat-pipe cooling than high angles of 
attack. A high angle of attack will heat a larger portion of the lower surface, making it less 
useful for radiating away heat transferred from the stagnation region. The heat must thus 
be moved to the upper surface which experiences very little heating. The required heat-pipe 
lengths are then much longer than for a correspondingly low angle of attack leading edge. 

The approximation for the design analyses presented here conservatively estimates the 
maximum temperatures by assuming no transfer of heat chordwise at the stagnation line. 
In Figure 4, the 2-D geometry is approximated as a 1-D problem. However, due to the 
sharp reduction hi heat flux at the stagnation line, three dimensions should be considered 
for a complete analysis. For sharp leading edges, heat will be transferred away from the 
stagnation line hi the chordwise direction parallel to the heat pipes, thus reducing the 
leading-edge temperatures and resulting hi a conservative approximation. For blunt leading 
edges, the chordwise heat flux reduction is much less, and the three dimensional effect is 
correspondingly less significant. 

A second conservative feature hi the developed design analysis (and the FEA) is the 
use of a constant applied heat flux on the outer surface with only radiation losses, hi actual 
aerodynamic heathig, the reduction hi the applied heat flux with rising surface temperature 
is much greater if convection to a hot surface is considered rather than assuming a constant 
heat flux with radiation losses alone. This effect is most pronounced as the heat-pipe 
spacing hicreases. For a very small heat-pipe spacing, the surface temperature, and thus 
heat flux, is relatively uniform. For a large heat-pipe spacing, the maximum temperature 
between heat pipes is much greater than directly over a heat pipe. If surface temperature 
dependent convective aerodynamic heathig is considered instead of a constant applied heat 
flux, the aerodynamic heathig applied to the surface will decrease significantly with the rise 
in surface temperature. This results hi a “damping” of the temperature rise. The design 
equations presented here can thus be used to evaluate feasibility assuming closely spaced 
heat pipes. Once it is determined that heat pipes are feasible with closely spaced heat pipes, 
a more detailed analysis should be utilized to determine optimum heat-pipe spachig. 

Concluding Remarks 

A set of closed form equations have been presented to quickly evaluate the feasibility 
of utilizing heat pipes to cool leading edges of hypersonic vehicles. The results from the 
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design equations were compared with results from a 3-D FEA for both a large and small 
radius leading edge. The results compared quite well and indicate that the equations can be 
used for a quick assessment of the feasibility of using heat pipes to cool a leading edge. If 
feasibility is indicated, a more detailed analysis should follow. 
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Appendix A: Spreadsheet Description 

A spreadsheet was used to determine the temperatures using the design equations. An 
example of the spreadsheet used in shown in Table 4 for test case 2. Rows 2-15 are the 
input variables, many of which were listed in Table 1. For the thennal conductivity and 
heat flux per unit area, the values given in column D have the units listed with the definition 
in column B. The values in column E have the units listed in column F. Rows 17-21 are 
the output values. The parameter, along with its units, is listed in column B and the 
calculated value is listed in column C. The equations used to calculate the values in column 
C are given in column D. 


Table 4: Spreadsheet Showing 1 

Leading-Edge Temperature Calculations 

1T1I 

B 

C 

D 

E 

F 

2 






3 

t wall, in. 



0.01 


4 

t coating, in. 



0.01 


5 

half hp spacing, in. 



0.05 


6 

L upper, in. 



36 


7 

L lower, in. 



18 


8 

k structure, t-t, Btu/hr-ft-°F 

12.8 

0.000296 

Btu/in-s-°F 

9 

k structure, plane, Btu/hr-ft-°F 

24.408 

0.000565 

Btu/in-s-°F 

10 

kwall, Btu/hr-ft-°F 

39.485 

0.000914 

Btu/in-s-°F 

11 

k coating, Btu/hr-ft- 

°F 

23.587 

0.000546 

Btu/in-s-°F 

12 

emittance 



0.8 


13 

q stag, Btu/ft2-s 


83 

0.576 

Btu/in2-s 

14 

qtot, integrated heat load, Btu/s 


12.09 


15 

heat-pipe width, in. 



0.6 
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17 

R 

242.41 

=E4/E11 + 
E4*E11) 

E2/E8 + E5 *(E2+E4)/ (E2*E9 + 

18 

AT stag, °F 

140 

BslfcHIM 



19 

Tsurf, °F 

2573 

=(E14*3600*144/(E12*0.000000001714*(F6+ 
E7))) A 0.25 - 460 

20 

AThp, °F 

43 

=((E14/((E6+E7)*(E15/(E15+2*E5)))))*(E2/E8 
+ E3/E10 + E4/E11) 

HI 

Thp, °F 

2530 




22 

Tmax, °F 

2670 

guniBi 
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Appendix B: Finite Element Analysis 

A three-dimensional, thermal finite element model was used to obtain the finite 
element solutions to compare with the design results. A schematic diagram of the leading 
edge is shown in Figure 9. The shaded region in the figure represents surfaces of the 3-D 
region that was modeled in the analysis. The Engineering Analysis Language (EAL) 
system was used to perfonn the finite element analysis. 16 The model was constructed 
using the executive control language hi EAL in a very general sense, in that the physical 
dimensions of the leading edge, the composite architecture, and the boundary conditions 
can be easily varied. 



Figure 9: Schematic drawing of a wing leading edge showing surfaces of the 3-D region 
modeled in the finite element analysis. 

Elements were clustered near the stagnation region in the chordwise direction where 
the heat flux and temperature gradients are the largest. Element clustering was used in the 
chordwise direction since the chordwise dimension is much larger than the spanwise or 
through-the-thickness dimensions. The transformation used to concentrate the elements hi 
the stagnation region hi the chordwise direction is a logarithmic clustering algorithm given 
by 


where 


y = B + 


X 


-1 

sinh 


— - 1 ]sinh(xB) 
Vz ) 


1 

B = — In 
2t 


l + (e T — l)(z / s) 
l + (e _T -l)(z/s) 


, 0 < x < 


oo 


( 8 ) 


(9) 


and where y is the transformed coordhiate in the chordwise direction along the leading 
edge, p the orighial equally spaced coordhiate, z the pomt about which the clustering 
occurs, and s the heat-pipe length. The stretching parameter x can be varied to space the 
pomts equally (small x) or to concentrate the points near z (large x). In the spanwise and 
through-the-thickness directions a linear grid was used. A typical finite element model 
used for the comparison is shown hi Figure 10. In the figure, the chordwise length is 
much smaller than the actual case, but the number of elements is representative of those 
used in the analysis. 
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Figure 10: Finite element model of section of leading edge modeled. 

The walls of the heat pipes were modeled by three-dimensional conduction elements 
and the internal vapor temperature was modeled by an isothermal surface on the hmer 
surface of the heat pipe to simulate an infinite thermal conductivity for the heat pipe. This 
was a nonconservative assumption, and detailed heat-pipe analyses are necessary to 
determine actual temperature drops along the heat pipe. However, it is a good assumption 
for high-temperature liquid metal heat pipes. 

The cross-sectional finite element grid used in the model is shown in Figure 1 1 along 
with the boundary conditions. Radiation exchange between the hot leading edge and a cool 
ambient is insignificantly different than if the ambient is assumed to be at absolute zero, and 
thus heat is radiated from the external surface to space at absolute zero. (This assumption 
is non-conservative, but has a negligible impact on the surface temperatures with high 
aerodynamic heating.) An insulated internal surface was considered since all interior 
surface temperatures are relatively uniform. The heated surface and the surfaces cooled by 
radiation are shown hi Figure 1 1 . The surface midway between the chordwise heat pipes 
and the surfaces through the center of the chordwise heat pipe are assumed to be thermally 
insulated as a result of symmetry. 

Several simplifying assumptions were made in the finite element analysis. Perhaps 
the largest uncertainty hi the thennal finite element analysis is the assumption of perfect 
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thermal contact between the carbon-carbon structure and the heat pipes. It is known that 
thermal contact resistance will result hi increased surface temperatures hi the stagnation 
region, but smce no value for the contact resistance is known, its quantitative effect on the 
temperatures is uncertain. It is anticipated that, upon heating, the thermal expansion of the 
heat pipe would mcrease the contact pressure between the heat pipe and structure, thereby 
reducing the thermal contact resistance. 
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Figure 11: Boundary conditions used in the finite element analysis. 
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